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The impact of an air-sea roughness parameterization over the ocean that 
more closely matches recent observations of air-sea exchange is examined in 
the NASA Goddard Earth Observing System, version 5 (GEOS-5) atmospheric 
general circulation model. Surface wind biases in the GEOS-5 AGCM are 
decreased by up to 1.2m/s. The new parameterization also has implications 
aloft as improvements extend into the stratosphere. Many other GCMs (both 
for operational weather forecasting and climate) use a similar class of param- 
eterization for their air-sea roughness scheme. We therefore expect that re- 
sults from GEOS-5 are relevant to other models as well. 
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1. Introduction 


The interaction between the ocean surface and the lowest levels of the atmosphere is a 
crucial component of any atmospheric GCAL The exchange of momentum, moisture, and 
sensible heat between the ocean and atmosphere occurs on spatial and temporal scales 
far finer than any GCM can directly simulate. Many models therefore rely on Monin- 
Obhukov Similarity Theory (MOST) to specify air-sea exchange as a function of bulk 
winds, temperature, and humidity. Early attempts at quantifying the exchange coeffi- 
cients underlying MOST were conducted under conditions far removed from that actually 
experienced in the ocean (e.g. Charnoek [1955]; Large and Pond [1981]). Nevertheless, 
generations of atmospheric models have relied on these earlier measurements for tuning 
their air-sea roughness scheme. For example, G EOS-5 currently implements Large and 
Pond [1981] for moderate and strong winds and Kondo [1975] for weak winds [Helfand 
and Schubert, 1995]. See Table 1 for a description of the schemes in a range of models. 

More recent in-situ observations have improved our understanding of air-sea exchange 
over deep ocean waters, especially over high wind regions like the Southern Ocean. In 
particular, recent field campaigns have measured turbulent exchange over the Southern 
Ocean, over the Gulf Stream, and over the North Atlantic in high wind speeds (e.g. Edson 
[2008], Edson, in preparation for Journal of Physical Oceanography, Yelland et al. [1998], 
Edson et al. [2007], and Banner et al. [1999]). These field campaigns have found that 
the Charnok parameter appears to increase with wind speed beyond lOrn/s, so that a 
parameterization based on Charnoek [1955] or Large and Pond [1981] underestimates the 
drag on surface winds (section 3c of Faira.ll et al. [2003]). Recent observations of air-sea 
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33 exchange imply that the current air-sea roughness scheme in GEOS-5 produces too little 
drag on surface winds in the range of wind speeds common in the Southern Ocean, 
k Accurate climatologies of surface winds over ocean regions were not available when the 
56 current Large and Pond [198.1]-basecl parameterization in the GEOS-5 model was created, 
sr but satellite-based climatologies of surface winds are now available (e.g. Chou et al. 
.a [2003]). These satellite based climatologies suggest that surface winds in the GEOS-5 
33 model are too strong over the Southern Ocean and off the coast of Asia in the North 

40 Pacific (Figure la-c). As surface winds over the Southern Ocean drive present and future 

41 oceanic uptake of C0 2 [Downes et ah, 2011; Matebr and Hirst, 1999], it is important to 
45 accurately simulate surface climate in this region. The GEOS-5 model is not alone in its 

43 poor representation of Southern Ocean surface wind; Barnes and Hartmann [2010] find 

44 that the latitude of the Southern Hemisphere jet maxima varies by over 5° in the Coupled 

45 Model Intercornparison Project (CMIP3) ensemble, and that such a bias has implications 

46 for the response of a GCM to doubled C0 2 - Regional models also have difficulty capturing 

47 mesoscale turbulent surface fluxes [Renfrew et ah, 2009]. 

48 This paper discusses efforts to reduce this bias in GEOS-5 by updating the air-sea 
43 roughness parameterization from Helfand and Schubert [1995]. Section 2 describes changes 
so made to the model and Section 3 presents results. As other atmospheric GCMs base their 
si air-sea roughness parameterization for momentum exchange on similarly old data, we 

52 expect that the reduction in model bias shown here might be common to other GCMs as 

53 well. 
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2. Change to Scheme 

We first describe the air-sea roughness scheme in G EOS-5 before discussing the changes 
made to increase the surface friction. G EOS-5 contains 72 vertical levels, with approx- 
imately 8 in the boundary layer fet al., 2008]. Like many atmospheric GCMs, G EOS-5 
uses MOST to describe momentum, heat, and moisture flux coefficients in terms of bulk 
quantities (e.g. zonal wind, specific humidity, and temperature) in the model. The wind 
stress vector at the surface can be expressed as 

[t x , t v ] = p a v $ C D A[u, w]; |r| = p a u* 2 , (l) 

where p a is the air density, Co the transfer coefficient for momentum, v s the difference be- 
tween the ocean and atmosphere surface wind speed, A[w, v] is the difference between the 
ocean and atmosphere surface wind vector, and u* the friction velocity. MOST computes 
u* (and Co) as a function of bulk parameters via the following equations: 

Co = k~ jdG/o(C~ : o)] ' , 

u* — Cp 2 v (2) 

Zq — + Ai T A$U * -\-A4U * 2 

54 where zq is the roughness length, k is the Von-Karman constant, Ai through An are tunable 

55 parameters used to match the air-sea roughness scheme to observations, and $mo( C 2 o) is 

se controlled by stability of the air column above. After an initial guess is made at Cp (in 

57 practice Co assuming neutral stability), Equation 2 is solved iteratively until a new value 

sa for Co has been reached consistent with the actual stability. 

53 Previously, the Ai through A 5 coefficients were chosen to interpolate between the recip- 

60 rocal relation of Kondo [1975] for weak winds and the piecewise linear relation of Large 

6 1 and Pond [1981] for moderate to large winds. The key change described by this paper 

62 is that the values for the Ai through ,4 5 coefficients are changed so that the roughness 
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63 length is increased for a given friction velocity. Neither the formulation for TuoftAo) nor 
6, the coefficients at low wind speeds is changed. For very strong winds (e.g. hurricanes), 
es roughness length no longer increases with wind speed [Molod and Partyka, 2011]. See 
66 Table 1 for the coefficients used. Runs with the old polynomial for : 0 are referred to as 
s? CONTROL, and runs with the new polynomial for are referred to as NEW. 
es Figure 2a shows Con mm as a function of 10m wind speed for the old and new coefficients 
es and in observations. The drag coefficient has been increased beyond the average suggested 
ro by the most recent observations but within the uncertainty. We chose the highest drag 
?! coefficient justified by the observations to achieve the maximum impact on the GEOS-5 

72 wind bias. Any further increase would distance us from the range of observational nncer- 

73 tainty. Note that the drag coefficient in Community Atmosphere Model (CAM/CCSM) 

74 (dashed red line) appears to be too small. CAM has not upgraded its scheme since version 

75 2.0 (Kiehl et al. [1998] versus section 4.11.2 of Neale and et al. [2010]). Figure 2b compares 

76 modeled output roughness length and friction velocity for the NEW and CONTROL runs. 

77 As expected, surface roughness dramatically increases with the new coefficients. Figure 

78 2b also includes curves of zq = achamok'U* 2 / 9 [Charnock, 1955] but with different values of 

79 the Charnok parameter achamok- Older measurements suggest values of achamok ~0.011 
so to achamok ~0.018 (see Section 3c of Fairall et al. [2003]). Newer observations (Edson 
si [2008] and Edson, manuscript in preparation for .Journal of Physical Oceanography) would 
»2 imply a higher a C hamok- 
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83 This change has been implemented in GEOS-5. Several GEOS-5 atmosphere-only sim- 
» illations with the old and new coefficients were performed to examine the impact of the 
as increased drag: 

86 1. 2x2.5 degree 30 year runs with interactive stratospheric chemistry, 

a? 2. 2x2.5 degree 12 year run without interactive stratospheric chemistry, 

8 a 3. 1x1.25 degree 25 year run without interactive stratospheric chemistry, 

as 4. a series of 1/4 degree 5-day forecasts. 

<w All simulations showed similar impact of the new' roughness parameterization, and we 
oi will focus here on results from the 30 year run with stratospheric chemistry. CONTROL 
92 and NEW differ only in the air-sea roughness scheme; all other models settings are fixed, 
os A Student- T two-tailed test is used to assess statistical significance. Each year is taken as 

94 one degree of freedom. Surface winds and surface stress from Version 2 of the Goddard 

95 Satellite-Based Surface Turbulent Fluxes (GSSTF ) Data [Chou et ah, 2003] are used to 

96 validate the model. We now address the impact of this change iri the air-sea roughness 
•)? parameterization on bulk quantities in the model. 

3. Results 

98 We now discuss how the change in friction influences the momentum budget in the 

99 model. The exchange coefficient for momentum increases over most oceanic regions, with 

100 the strongest increase over the Southern Ocean (Figure 3b). Biases in surface wind are 

101 reduced across the ocean regions in response to the altered surface roughness coefficients 
M 2 (Figure lc-d). Winds over the Southern Ocean decrease by over lrn/s, but winds are 
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.os reduced over most ocean covered regions. Surface sea level pressure biases are also reduced 
io« (not shown), consistent with the wind speed improvement. 

,03 Figure 3c-f shows zonal surface stress on the ocean. Changes in surface stress are 
m smaller than changes in either Co or wind speed, as might be expected from Equation 
is? 1. Namely, the decrease in wind speed and increase in Co largely balance each other, 
io 8 so that their product is nearly constant. Nevertheless, the changes are significant in the 
us Southern Hemisphere, whereby surface stress on the Southern Ocean is increased while 
no surface stress further equatorward is decreased. The change is particularly strong in the 
m Indian Ocean/ Australia region. Biases in the control run are partially ameliorated. Runs 
no in which the atmosphere is coupled to a full ocean are planned in order to understand the 
us potential impact on the ocean circulation. 

lM These changes in surface stress imply anomalous eddy momentum flux convergence aloft, 
us as vertically averaged must balance surface friction for a steady state surface jet ( Held 
ns [1975] and section 12.1 of Vallis [2006]). Figure 4 shows that poleward momentum flux 
nr is increased throughout the upper troposphere, as implied by the dipole of surface stress, 
us Eddies are fluxing more momentum poleward in order to counteract the weakening of the 
us surface jet. Associated with this change in momentum flux are statistically significant 

120 improvements in extratropical forecasting skill (not shown). 

4. Conclusions 

1 2 1 The old air-sea roughness scheme in GEOS-5 is based on 30-year old observational 
,22 data, but newer data suggests seas are rougher. Associated with the old parameterization 
123 are overly strong surface winds. By incorporating more recent observations of air-sea 
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S 3 » exchange into the model's air-sea roughness scheme, we have improved the surface climate 
’25 in the G EOS -5 AGCM, Preliminary results indicate that tire improvement is present at 
126 resolutions up to 1/4 degree. 

is? Modifying the air-sea roughness parameterization leads to statistically significant 
ns changes in cloud distribution, heat flux, stratospheric ozone, and planetary wave driving 
l29 of the stratosphere. Presentation of these changes, a discussion of the surface moisture 
uo and sensible heat budgets, and further diagnostics on the tropospheric momentum budget, 

1 3 1 will be reported in detail in a future paper. The microphysics scheme in all runs con- 

132 sidered does not include interactive aerosals; preliminary results indicate that including 

133 interactive aerosals along with this change in surface roughness leads to large changes in 

134 sea salt aerosal concentration and subsequent cloud formation. 

135 Other atmospheric GCMs appear to use a similar scheme to parametrize the exchange 
ne of momentum, heat, and moisture with the ocean. We expect that biases in these other 

137 models might be reduced if these models were retuned to more closely match available 

138 observations. 
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wind speed, ann. avg. 


(a} Control (b) Obs 



C.I.. 0 7m/s 


Figure I. (a) surface wind speed in the control run, (b) surface wind speed in observations, 
(c) control minus the observations, and (d) the new run minus the control. For (a) and (b), the 
contour interval is 2m/s and the color scale is on the top left. For (c) and (d), the contour interval 
is 0.7 m/s. For (c), the color scale is on the left. For (d), regions with anomalies whose statistical 
significance exceeds 95% are in color. The zero contour is omitted and negative contours are 
d fished. 
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Drag Coefficients from Observations, COARE, and GEOS-5 MOST relationship, GEOS-5 



o 5 10 15 20 25 o 0.2 0.4 0.6 0.8 1 


U ion (m/s) u* (m/s) 

Figure 2. (a) Neutral drag coefficient for momentum exchange at the ocean surface ( Coniom ) 

as a function of wind speed at 10m in observations [Banner et al., 1999: Yelland et ah, 1998) 
and in models. COARE3.0, COARE4.0, ECMWF wave model (i.e. not the uncoupled atmo- 
spheric model as in Table 1), and binned data are based on Edson [2008] and Edson. personal 
communication. Error bars for binned data denote 1 standard deviation. Model results are from 
CAM2.0-CAM5 [Kiehl et ah, 1998], and the original and new curves from GEOS-5. (b) Relation- 
ship between friction velocity (u*J and roughness length(zo) over all ocean gridpoints averaged 
over one day of GEOS-5 model output. Isolines of r 0 = achamok'U * 2 j g [Charnock, 1955] but 
with different values of the Charnok parameter atchamok are included for comparison. 
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jj) Ne w- Control, ‘M mem e x change 


4 , mom exchan 


C.I.: 0.001 kgm'V 2 

• d ) 0^ xtau eastward surface stress 


(c) Control, xtau, eastward surface stress 


(f) New- Control, xtau, eastward surface stress 


(e) Control-Obs, xtau, eastward surface stress 


Figure 3. (a,b) Cm, drag coefficient for momentum exchange at the surface (C' m =Co*surfa.ce 


in the control run and in the new run minus the control. Contour interval is 


5 • 10 3 kgm 2 s 2 for (a) and 10 3 kgm 2 s 2 for (b). (c-f ) Eastward surface stress at the surface 


in the control run (c), observations (d), control-observations (e), and new-control (f). Contour 


interval is 5 • 10 2 Nm 2 for (c) and (d) and 10 2 Nm 2 for (e) and (f). Regions with anomalies 


whose statistical significance exceeds 95% are in color in (b) and (f). The zero contour is omitted 


and negative contours are dashed 
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(a) Control 


c) Control-Obe 


40 

latitude 


40 

latitude 


Figure 4. Momentum flux (■ u'v ') latitude- height cross section. Contour interval is 12.5m 2 s -2 
for (a) and (b), 2 m 2 s -2 for (c), and lm 2 s~ 2 for (d). The zero contour is omitted and negative 
contours are dashed. Regions with anomalies whose statistical significance exceeds 95% (95%) 
are in light(dark) blue in (d). 


D E A F T 


August 15, 2011, 6:04pm 


DRAFT 






H 

tL, 

< 

BS 

Q 


(a) Momentum exchange over oceans in different GCMs. 


vl deling Group 

source 

description 

CAM 2 

Kiehl et al. fl 998] 

based on Large et al. [1994] 

CAM5 

p.181 of Neale and et al. [2010] 

Large et al. j 1 994] 

GEOS- 5 (old) 

Helfand and Schubert [1995] 

Large and Pond [1981] 

AM2.0 

model development team [2004] 

Beljaars [1995], achamok — 0.018 

Y\TIF(MM5) 

page 72 of Skamaroek et al. [2008] 

based on Charnock [1955] 

uncoupled ECMWF 

page 38 of for Medium-Range Weather Forecasts [2010] 

based on Charnock [1955], cxchamok — 0.018 


ib) Coefficients for MOST Scheme equation relating u* to ; 0 (t 0 = 4- A-j -f- A 3 u * +A 4 u * 2 -Kl 5 u< 3 ) 





A 2 

-4;j 

A 4 

-4| 

control 

M* <0.0632456 

0.20303255 - 5 

0 

0 

0 

0 


fl.Q632456<u*<0. 381844 

-0.402451E-08 

0.239597E-04 

0.117484E-03 

0.191918E-03 

0.395649E-04 


0.38 1844 <u* 

-0.237910E-04 

0.228221E-03 

-0.860810E-03 

0.176543E-02 

0.784260E-04 

new 

u' : G;0. 0632456 

0.20303255 - 5 

0 

0 

0 

0 


0.0632456<u* 

-1.1 0245 IE-08 

0.1593E-04 

0.1E-03 

2.918E-03 

0.G95649E-04 


Table 


The Community Atmosphere Model (CAM) has not upgraded its scheme since version 2.1) (Kiehl et al. [1998] 
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Ph 

PS 
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versus section 4.11,2 of Neale and et al. [2010]). 
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